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A spatially extended Lotka-Volterra system of two competing species 
in the presence of two correlated noise sources is analyzed: (i) an external 
multiplicative time correlated noise, which mimics the interaction between 
the system and the environment; (ii) a dichotomous stochastic process, 
whose jump rate is a periodic function, which represents the interaction 
parameter between the species. The moment equations for the species den- 
sities are derived in Gaussian approximation, using a mean field approach. 
Within this formalism we study the effect of the external time correlated 
noise on the ecosystem dynamics. We find that the time behavior of the 1** 
order moments are independent on the multiplicative noise source. How- 
ever the behavior of the 2"'' order moments is strongly affected both by 
the intensity and the correlation time of the multiplicative noise. Finally 
we compare our results with those obtained studying the system dynamics 
by a coupled map lattice model. 
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1. Introduction 

Real ecosystems are influenced by the presence of continuous external 
fluctuations connected to the random variation of environmental parame- 
ters such as temperature and natural resources, which affects the system 
dynamics by a multiplicative nonlinear interaction pj . The spatio-temporal 
behaviour and the formation of spatial patterns became recently an impor- 
tant topic in hydrodynamics systems, nonlinear optics, oscillatory chemical 
reactions and in theoretical ecology [2] -[SI- During the last three decades a 
new approach, based on the use of moments, has been exploited to describe 
the behaviour of spatially extended system in population dynamics, quan- 
tum systems in the context of nonlinear Schrodinger equations, and kinetic 
models of polymer dynamics [5]-[7]. In this paper, by using the formalism of 
the moments, we study the spatio-temporal behaviour of a two-dimensional 
system formed by two competing species subject to random fluctuations. 
The system is described by generalized Lotka-Volterra equations in the pres- 
ence of two noise sources: (i) a multiplicative time correlated noise with 
correlation times t'^, modeled as an Ornstein-Uhlenbeck process [8], which 
takes the environment fluctuations acting on the species into account; (ii) a 
noisy interaction parameter which is a stochastic process, whose dynamics 
is given by a periodic function in the presence of a correlated dichotomous 
noise, with correlation time r^. We define a two-dimensional spatial domain 
considering in each site a system of two Lotka-Volterra equations coupled 
by an interaction term P]. Afterwards, using a mean field approach, we 
study the dynamics of the system by the moment equations, within the 
Gaussian approximation [10\ lllj . getting the time behavior of the 1** and 
2nd Qpjjgj. moments of the species concentrations. Finally we compare our 
results with those obtained within the formalism of the coupled map lattice 
(CML) model [l2]. 

2. The model 

Our system is described by a time evolution model of Lotka-Volterra 
equations, within the Ito scheme [I3]-[T5], with diffusive terms in a spatial 
lattice with sites 

7 

yi,j = mA'^ - yi,j - P^hi) + VhiV^Clj + DY^{y^- Vij), (2) 

7 

where Xij and Uij denote respectively the densities of species x and species 
y in the lattice site fJ- is the growth rate, D is the diffusion constant, 

and ^2-^ indicates the sum over all the sites except the pair Here 
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f5 is the interaction parameter. C (^) = ^lU) ^-re statistically indepen- 
dent colored noises, i. e. exponentially correlated processes given by the 
Ornstein-Uhlenbeck process [8] 

i = -^C + ^«'«) = (3) 

and ^'(t) (/ = are Gaussian white noises within the Ito scheme with 
zero mean and correlation function {i\t)£,'^{t')) = 2ad{t - t')Sim- The 
correlation function of the processes of Eq.(l3|) is 

(C'(i)r (t')) = ^e-l*-*'l^^(5z^ (4) 
and gives 25 {t — t')5im in the limit — > 0. 

2.1. The interaction parameter 
The value of the interaction parameter (3 is crucial for the dynamical 
regime of the ecosystem investigated. In fact, for /3 < 1 both species survive 
and a coexistence regime takes place, while for /3 > 1 one of the two species 
extinguishes after a certain time and exclusion occurs. These two regimes 
correspond to stable states of the Lotka-Volterra's deterministic model |13j . 
Moreover periodical and random driving forces connected with environmen- 
tal and climatic variables, such as the temperature, modify the dynamics 
of the ecosystem, affecting both directly the species densities and the in- 
teraction parameter. This causes the system dynamics to change between 
coexistence (/3 < 1) and exclusion {(3 > 1) regimes. To describe this dy- 
namical behavior we consider as interaction parameter P{t) a dichotomous 
stochastic process, whose jump rate is a periodic function 

lit) =0, At< Ta 

(5) 

7(*) = 70 (l + A \ cos ujt\) , At > . 

Here At is the time interval between two consecutive switches, and is 
the delay between two jumps, that is the time interval after a switch, before 
another jump can occur. In Eq. ([SJ, A and uj = {2'k)/T are respectively 
the amplitude and the angular frequency of the periodic term, and 70 is the 
jump rate in the absence of periodic term. Setting (3down = 0.94 < 1 and 
(3up = 1-04 > 1, the dichotomous noise causes f5{t) to jump between two 
values, Pdown and Pup- The value = 43.5 for the delay corresponds to 
a competition regime with /3 switching quasi-periodically from coexistence 
to exclusion regimes (see Fig. [1]). This synchronization phenomenon 
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Fig. 1. Time evolution of the interaction parameter f3{t), driven by time correlated 
dichotomous noise with correlation time Td = 43.5. P(t) jumps between (3down = 
0.94 and Pup = 1-04. The initial value is /3(0) = Pup = 1-04. The values of the 
other parameters are: A = 9.0, lo/{2tt) — 10~^, 70 = 2 • 10^^. 

is due to the choice of the value, which stabilizes the jumps in such a 
way they happen for high values of the jump rate, that is for values around 
the maximum of the function 7(t). This causes a quasi-periodical time 
behaviour of the species concentrations x and y, which can be considered 
as a signature of the stochastic resonance phenomenon [H] . 

3. Mean field model 

In this section we derive the moment equations for our system. Assuming 
N ^ 00, we write Eqs. ([T|) and ([2]) in mean field form 



where < x > and < y > are average values on the spatial lattice considered, 
that is the ensemble average in the thermodynamics limit. We set fxix, y) = 
fix{l - X - (5y), gx{x) = x, fy{x,y) = - y - (3x), gy{y) = y. By site 
averaging Eqs. ([6]) and we obtain 



By expanding the functions fx{x,y), gx{x), fy{x,y), gy{y) around the 1** 
order moments < x > and < y >, we get an infinite set of simultaneous 
ordinary differential equations for all the moments. This equation set is 
truncated by using a Gaussian approximation, which causes the cumulants 
above the 2"'^ order to vanish. Therefore we obtain 



X = fx{x,y) + y/axgx{x)C + D{< X > -x) 
y = fy{x,y) + ^/aygy{y)C + D{< y > -y), 



(6) 
(7) 




(8) 
(9) 



< X > = < X > (1— < X > —j3 < y >) — )u(/3/iii + /i2o) 



(10) 
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<y> = IX < y > {I- < y > < X >) - /u(/3/iii + /io2) (H) 
fi2Q = 2/X/X20 - 2D^2Q - 2/U/3 < y > iJ,2o - 2n < X > + 2//2o) 

+ 2a^{< X >2 +/i2o)(l - e"^) (12) 
/io2 = 2/X/X02 - 2I);Uo2 - 2(U/3 < X > IJ.02 - 2fi< y > + 2/io2) 

+ 2ay{< y >2 +;uo2)(l - e"^) (13) 
/ill = 2/z/xii - 2DIJ.U- <x> [2fiiJ,u + /x/3(/iii + /io2)] 

- < y > [2/i/iii + fi/3{nu + /i2o)] , (14) 

where ^20 : A'02, /^ii are the 2""^ order central moments defined on the lattice 

fi2o = <x^ > - <x >^ fjLQ2 =<y^ > - <y >^ (15) 
flu = < xy > — < X >< y > . (16) 

The dynamics of the two species is analyzed through the time behavior of 
1** and 2"'^ order moments according to Eqs. (110p - (jl4p . for different val- 
ues both of intensities a {= = (Jy) and correlation time r'^ (=7"x 
of the multiplicative colored noise. The results are reported in Figs. [2]| 
The results shown in Fig. [51 obtained with a low value of the correlation 
time t'^ reproduce the same time behavior obtained for uncorrelated white 
noise sources in Ref. [11] . The value of the time delay is = 43.5, which 
determines a quasi-periodic switching between the coexistence and exclu- 
sion regimes. The values of the other parameters are: fi = 2, D = 0.05. 
The initial conditions are: C(0) = C^(0) = 0, < x{0) > = < y{0) > = 
0.1, /U2o(0) = /Uo2(0) = /Uii(O) = 0. These initial values for the moments 
correspond to uniformly distributed species on the considered lattice. In 
Figs. [2^, [3^, we note that the 1*** order moments undergo correlated os- 
cillations around 0.5. This behaviour is independent both on the intensity 
and the correlation time of the multiplicative noise. On the other hand 
the behavior of the 2"*^ order moments depends strongly both on the in- 
tensity and the correlation time of the external multiplicative noise. In the 
absence of noise fj,20, /^02) l^ii maintain their initial values. For very low 
levels of multiplicative noise {a = 10~^^) quasi-periodical oscillations ap- 
pear with the same frequency of the interaction parameter /3(i), because 
the noise breaks the symmetry of the dynamical behavior of the 2"^^ order 
moments (see Figs. [2)3, [HJa, HJa). About the variances of the two species, 
the time series of H20 and /X02, which coincide all the time, show an os- 
cillating behaviour characterized by small (close to zero) and large values. 
However the negative values of the correlation //n indicate that the two 
species distributions are anticorrelated. In particular, we find a time behav- 
ior characterized by anticorrelated oscillations, whose amplitude increases 
with the multiplicative noise intensity and it is reduced as the correlation 
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Fig. 2. Time evolution of the P* and 2'"^ order moments for = 1. The time series 
of (a) < x{t) > and < y{t) >, and (b) ^2,0 and ^0,2 respectively, are completely 
overlapped. The values of the multiplicative noise intensity are: ct = 0, 10~^^, 10^^, 
from top to bottom. Here — 43.5, fJ. = 2, and D = 0.05. The initial conditions 
are: C"(0) = C^O) = 0, < x(0) > - < y(0) > = 0.1, ^^20(0) - ^02(0) - A^nW = 0. 
The values of the other parameters are the same of Fig. [TJ 

indicates that the spatial distribution in the lattice will be characterized by 
zones with a maximum of concentration of species x and a minimum of con- 
centration of species y and viceversa. The two species will be distributed 
therefore in non-overlapping spatial patterns. This physical picture is in 
agreement with previous results obtained with a different model ^X5\- A 
further increase of the multiplicative noise intensity {a = 10~^) causes an 
enhancement of the oscillation amplitude both in fj,2o, /X02 and fin. This 
gives information on the probability density of both species, whose width 
and mean value undergo the same oscillating behavior. The anticorrelated 
behavior is enhanced by increasing the noise intensity value. We note that 
the amplitude of the oscillations is of the same order of magnitude of the 
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Fig. 3. Time evolution of the 1** and 2 order moments for t'^ = 10. Here 
(a) < x{t) > and < y{t) >, and (b) ^2,0 and ^0,2 respectively, are overlapped. 
The values of the multiplicative noise intensity are: a = 10^^^, 10^^, from top to 
bottom. The values of the other parameters and the initial conditions are the same 

of Fig. m 



noise intensity a, that is the amplitude of the oscillations is enhanced as the 
noise intensity increases. The periodicity of these noise-induced oscillations 
shown in Figs. [2]|3]is the same of the interaction parameter (see Fig. [1]). 
Moreover the right-hand side of Figs. [2]|3] shows that the multiplicative col- 
ored noise affects the time evolution of the 2"^^ order moments introducing a 
delay: the amplitude of the oscillations reaches its highest value after a time 
interval whose length increases as becomes bigger. In fact by comparing 
the first figures in Figs. [2t)-|D3 with a = 10^^^, for example, the maximum 
is reached approximately at t ~ 110 for = 1, and at t ~ 330 for = 100. 

4. Coupled Map Lattice Model 

Our results, obtained by using moment equations in Gaussian approxi- 
mation, can be checked studying the dynamics of the system by a different 
approach, namely the CML model 

^Sr^ = M-g(l-4?-/3^"^yS) + ^/^^gC•^^+^E(-^"^-4?) (17) 

7 

which represents a discrete version of the Lotka-Volterra equations. 
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Fig. 4. Time evolution of the 1*** and 2"'' order moments for r'^ = 100. Here 
(a) < x{t) > and < y{t) >, and (b) fj,2,o and /Lto,2 respectively, are overlapped. 
The values of the multiplicative noise intensity are: a — 10^^^, 10^*^, from top to 
bottom. The values of the other parameters and the initial conditions are the same 

of Fig. m 



■"^ and ylj denote respectively the densities of preys x and y 



Here x 

in the site at the time step n, fj, is the growth rate and D is the 

diffusion constant. The interaction parameter /^^"^ corresponds to the value 



of P{t) taken at the time step n, according to Eq. Q^Y^ and Cf} 

are independent Gaussian colored noise sources above defined in Eq. ([S]). 
Finally indicates the sum over the four nearest neighbors. To evaluate 
the 1** and 2"'^ order moments we define on the lattice, at the time step n. 



the mean values < x >^A. , 

CML ' 



< z> 
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(19) 



the standard deviations var^\ vary^^ 
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Fig. 5. (a) Mean values: < x >cml^ < V >cmlj ^i^^ (b) variances: var^, vary 
and covxy of the two species, as a function of time for = 1. The values of the 
multiplicative noise intensity are: a = 0, 10~^^, 10^^, from top to bottom. The 
initial values of the species concentrations are xf'j = yfj =0.1 for all sites 
The values of the other parameters are the same of Fig. ^ 



and the covariance covmj of the two species 

covf^ = — ^ ^ ^ , (21) 

where = 100 x 100 is the number of lattice sites. We note that < x >(^j^jj;^, 

< y >CML varx, vary, covxy, corresponding respectively to < x >, 

< y >, and iJ,2o, fJ'02, f^ii obtained within the mean field approach, are the 
I'** and 2"^^ order moments defined within the scheme of the CML model. 
The time behavior of these quantities, for three different values of the mul- 
tiplicative noise intensity, and for three different values of the correlation 
time is reported in Figs. [SHTl Comparing these results with those shown 
in Figs. [2]S1 we note that the time behaviour obtained for the 1*** and 2"^^ 
order moments by the CML model is in a good qualitative agreement with 
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Fig. 6. (a) Mean values: < x >cml^ < V ^cm-li ^-iid (b) variances: var^^ vary 
and covxy of the two species, as a function of time for t'^ — 10. The values of 
the multiplicative noise intensity are: a = 10^^^, 10^^, from top to bottom. The 
values of the other parameters and the initial conditions are the same of Fig. [5l 




Fig. 7. (a) Mean values: < x >cmli ^ V ^cml^ ^'^'^ (b) variances: var^, vary 
and coVxy of the two species, as a function of time for r'^ = 100. The values of 
the multiplicative noise intensity are: a = 10~^^,10~^, from top to bottom. The 
values of the other parameters and the initial conditions are the same of Fig. [5l 
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those found using the mean field approach. In particular we note the role 
played by the multiplicative noise: higher noise intensity causes the 2"'^ or- 
der moments to increase (see in Figs. [5117] the enhancement of the oscillation 
maxima as the multiplicative noise intensity increases). This indicates both 
a spread of the two species concentrations and a spatial anticorrelation be- 
tween them. These results agree with those found within the formalism of 
the moment equations. We observe that the time behaviour of varx, vary, 
covxy is characterized by the presence of oscillations whose maximum ampli- 
tude is reached after a time delay. This peculiarity is the same of that found 
within the mean field approach (see paragraph 3). The discrepancies in the 
oscillation intensities are due to: (i) the Gaussian approximation in the mo- 
ment formalism; (ii) the fact that the species interaction in the CML model 
is restricted to the nearest neighbors; (iii) the different stationary values of 
the species densities in the considered models. Specifically, using for j3 an av- 
erage value, (5aver = (Pup + Pdown) f^, We get: Xst = Vst = Paver) ^ 0.5 

for the mean field model, and x^^ = Vst = (1 ~ + Paver) — 0.25 

for the CML model. 



5. Conclusions 



By using the moment formalism in Gaussian approximation, we describe 
the time behavior of two competing species inside a two-dimensional spa- 
tial domain in the presence both of a multiplicative colored noise and a 
dichotomous noise. We find that the 1*** order moments of the two species 
densities show correlated oscillations, whose amplitude is independent on 
the multiplicative noise. However, the behavior of the 2"'^ order central 
moments depend strongly both on the intensity a and the correlation time 
t'^ of the multiplicative noise. In particular, the behavior of the 2"'^ order 
mixed moment /Un indicates that higher values of the multiplicative noise 
intensity push the two species towards an anticorrelated regime character- 
ized by oscillations whose maximum amplitude is reached after a delay time: 
this delayed behaviour depends on the correlation time t^. We find a good 
qualitative agreement between these results, obtained within the mean field 
approach, and those found by the CML model. In view of some applica- 
tions of our model, to describe and to predict the behaviour of biological 
species, we note that in real ecosystems the fluctuations are characterized 
by a cut-off. Therefore experimental data [16], whose dynamics is strongly 
affected by noisy perturbations and stochastic environmental variables, can 
be better modeled using sources of colored noise. 
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